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Correlation functions in lattice formulations of quantum gravity* 
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We compare different models of a quantum theory of four-dimensional lattice gravity based on Regge's original 
proposal. From Monte Carlo simulations we calculate two-point functions between geometrical quantities and 
estimate the masses of the corresponding interaction particles. 



1. INTRODUCTION 

There are two related schools attempting at 
a lattice theory of quantum gravity, dynamical 
triangulation and Regge quantum gravity. We 
concentrate on the latter approach and refer to 
for excellent reviews of the former frame- 
work. In particular we investigate different for- 
mulations, all of them based on the original work 
by Regge ||, namely conventional Regge grav- 
ity fl'll' a group theoretical approach ||, Z 2 - 
link Regge-theory |Q , and Regge gravity coupled 
to SU(2)-gauge theory ||. Each model exposes 
phase transitions at some critical gravitational 
couplings separating small and large curvature 
phases and allows to look for a continuum limit. 
|9HlCt|. A candidate realistic quantum theory of 
gravity should reproduce the expected long range 
interaction behavior observed in nature. Thus 
we calculate two-point funtions to probe the exis- 
tence of massless quanta of the gravitational field. 

2. LATTICE QUANTUM GRAVITY 

Any smooth d-manifold can be approximated 
by appropriately glueing together pieces of flat 
space, called d-simplices, ending up with a sim- 
plicial lattice. We take the edge lengths as degrees 
of freedom and leave the triangulation of the lat- 
tice fixed. Adopting the Euclidean path integral 
we can write down the partition function 



D[q, U]e 



-I{q,U) 



1) 



with an action / that covers all of the models 
introduced in more detail below. The functional 
integration extends over the squared edge lengths 
q and (if switched on) over the non-Abelian gauge 
fields U. 

One of the problems with (|l|) is the ambiguity 
in performing the link-lengths integration. Com- 
monly the measure is written as 



D[q]=Y[dq l qr 1 Hq) 



(2) 
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with T a function of the squared edge lengths q 
being equal to one if the Euclidean triangle in- 
equalities are fulfilled and zero otherwise. The 
question remains whether such a local measure 
is sufficient and how the power a might be cho- 
sen. A recent calculation for two-dimensional 
pure gravity pointed out that simple measures 
like (|^) do not respect the infinite volume of 
the diffeomorphism group. Only by appropri- 
ate gauge fixing and including the corresponding 
Faddeev-Popov term the correct continuum limit 
can be obtained pj| . However, the generaliza- 
tion of this procedure to higher dimensions and 
its numerical implementation are technically de- 
manding. 

Working in Euclidean space, i.e. with posi- 
tive definite metric, the conformal mode ren- 
ders the continuum Einstein-Hilbert action un- 
bounded from below. This unpleasant feature 
persists in the discretized Regge-Einstein action 
but need not necessarily lead to an ill-defined 
path integral Q|. Indeed, numerical simulations 
reveal the existence of a well-defined phase with 
finite expectation values within a certain range of 
the bare Planck mass, {m~) 2 < m 2 P < (m£) 2 
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On a triangulated lattice an action is given by 
I(q,U) = -2m%J2RMt,St) + Xj2 Vs 







^£V t Re[Tr(l-£7 t )] 



(3) 



The first sum runs over all triangle areas A t (q) 
and corresponding deficit angles St (q) yielding the 
curvature Rt. A cosmological term consisting of 
the cosmological constant A times the sum over 
the volumes V s (q) of all four-simplices s follows. 
Finally, an additional non-Abelian gauge action 
composed of the inverse gauge coupling (3, the 
weight factors W t {q) and the ordered product U t 
of SU(2) matrices around the triangle t is ap- 
pended. The weights 



Wt — const x 



A\ 



(4) 



with a four-volume Vt assigned to every trian- 
gle, describe the coupling of gravity to the gauge 
fields. 

Now we are ready to define the four differ- 
ent models for subsequent numerical treatment. 
Monte Carlo simulations have been performed 
on regularly triangulated hypercubic lattices with 
toroidal topology and 4 3 x 8 vertices. The gravi- 
tational couplings mp were chosen close to (m~) 2 
where there is a certain chance for a continuous 
phase transition ^] . 

2.1. Conventional Regge gravity 

Here we employ the Regge-Einstein action, 
putting Rt = Af5t with the gravitational cou- 
plings mp = {-0.0775, -0.0785, -0.0795k and 
include a cosmological term with A = 1 . No 
additional gauge fields are present in the action, 
(3 — 0, and the gravitational measure is chosen to 
be uniform, a = 1. 

2.2. Group theoretical approach 

Constructing the dual of a simplicial lat- 
tice, Poincare transformations can be assigned 
to its links to yield an action in which the 
sin of the deficit angle enters, Rt = A t sin St 
Q. Again we use a cosmological constant, 
A = 1, vanishing inverse gauge coupling, (3 = 0, 
the uniform measure, a = 1, and vary 
m% = {-0.055, -0.0555, -0.056}. 



2.3. Z2-link Regge-gravity 

This model is defined by restricting the squared 
link lengths to take on only two possible values 



1 + eai , oi E Z 2 , e < e J( 



6 M 



+ ■ 



(5) 



Then the quantum-gravity path-integral can be 
rewritten as the partition function of a spin sys- 
tem with somewhat complicated, yet local spin 
interactions 

2.4. Gauge fields coupled to Regge gravity 

In the system of SU(2) gauge fields cou- 
pled to quantum gravity we set Rt = A t St, 
A = and use a scale invariant measure, a = 
The following pairs of gravitational 

and gauge coupling, respectively, are considered: 
imp, (3) = {(-0.0025, 1.6), (-0.005, 1.0)}. 

3. TWO-POINT FUNCTIONS 

In order to examine the physical relevance of 
the well-defined phase mentioned above, we com- 
pute correlation functions of certain geometrical 
quantities. The gravitational weak-field propaga- 
tor can be cast into a spin-two and a spin-zero 
contribution |12|. The volume correlations are 
sensitive to the scalar part 

G v (d) = (j2v s J2 v »')c ( 6 ) 

sDv s'Dv d 

and the curvature correlations are due to the pres- 
ence of spin-2 particles 

tDv t'Dv d 

The local operators in (^|) and @) should be mea- 
sured at two vertices vq and Vd separated by the 
geodesic distance. We take the distance d to be 
equal to the index distance along the main axes 
of the skeleton. This seems a reasonable approx- 
imation in the well-defined phase with its small 
average curvature. In general one expects for @ 
and (0) at large distances the functional form 

p — rad 

A power law with a = 2 and a vanishing effective 
mass m = would hint at Newtonian gravity with 
massless gravitons. 



Rt / , Rf)c 



(7) 
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(a) Conventional Regge gravity: 
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(b) Group theoretical approach: 
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(c) Gauge fields coupled to gravity: O . . . {nip 
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Figure 1. Volume (left plots) and curvature (middle with magnification in the right plots) correlation 
functions for (a) conventional Regge gravity, (b) the group theoretical approach, and (c) the system of 
non-Abelian gauge fields coupled to Regge gravity. Error bars not explicitely drawn are in the size of the 
symbols. The curves correspond to fits with the function (8). 
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Table 1 

Effective masses m of Gr for several gravitational couplings mp 

(a) Conventional Regge gravity (b) Group theoretical approach (c) Gauge fields 

mp -0.0775 -0.0785 -0.0795 -0.055 -0.0555 -0.056 -0.0025 

m 2.9 3.1 3.0 1.1 0.9 0.5 2.1 



4. RESULTS 

Fig. 1 displays our Monte Carlo data of the 
two-point functions (|^) and (0). The volume cor- 
relations have already been studied around the 
transition at positive coupling where fits to an 
exponential decay have been obtained Jl3| . Such 
a fit procedure seems to be more difficult at neg- 
ative couplings. 

The curvature correlations are more suitable 
for a fit with (@). In order to test whether they 
obey a power law we fixed a = 2 and fitted the 
effective masses m. We took only distances d > 2 
into account. Gn(d = 1) is presumably plagued 
by lattice artifacts due to contact terms We 
are anyhow interested in the large distance be- 
havior. 

Table 1 contains the obtained mass parame- 
ters m. For (a) conventional Regge gravity m 
stays rather constant towards the critical coupling 
whereas in (b) the group theoretical approach the 
mass decreases for nip — ► (m~) 2 . In the case 
of SU(2) fields on the fluctuating lattice we get 
only one reasonable fit again indicating a nonzero 
mass. For all fits the uncertainties in the mass 
parameters are large. The Z2-link approximation 
of Regge gravity seems to behave differently. For 
this model the curvature correlations are compat- 
ible with a power law. The results of a more ex- 
tended study of the spin approach are reported 
in @. 

5. SUMMARY 

We computed two-point functions close to the 
critical bare Planck mass in the negative gravita- 
tional coupling regime. Altogether they exhibit 
a very similar behavior for the considered mod- 
els. Except for Z2-link Regge-gravity we found 
no convincing evidence for long range correlations 
corresponding to massless spin-zero or spin-two 



excitations. One reason might be, that contrary 
to the models described here, the ^-link theory 
is well-defined for all couplings mp and compu- 
tationally much less demanding, which allows to 
perform simulations exactly at the critical point. 
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